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Abstract. We calculate survival probability of a special state which couples 
randomly to a regular or chaotic environment. The environment is modelled by a 
suitably chosen random matrix ensemble. The exact results exhibit non-perturbative 
features as revival of probability and non-crgodicity. The role of background 
complexity and of coupling complexity is discussed as well. 
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1. Introduction 

The stability of a special prepared quantum state weakly coupled to a continuum is 
a subject of considerable interest in quantum information theory pQ, nuclear physics, 
mesoscopic, quantum chaos (see [H [3] and references therein). Giant resonances are 
collective excitations of the nucleons, which are approximate but not exact eigenstates 
of the complicated many-body Hamitonian. They are an example of a special state 
coupled weakly to a continuum picked from nuclear physics. Constructed superscar 
states in chaotic quantum billiards, considered recently [I] fall into the same class. The 
special state might also be implemented mesoscopically, for instance by an electronic 
state on a small conducting island (quantum dot), which is weakly coupled to one (or 
in a non-equilibrium situation to two) continua of the leads [5]. Given a particle in a 
quantum state ^(O)) at time t = the likelyhood to find the particle after some time 
t in the same state is measured by the survival probability 



The simplest approximation of P(t) is found perturbatively by F(ermi's) G(olden) R(ule) 
P(t) = e~ rt , where T is the inverse decay time, related to the mean coupling strength 
of the special state to the background. Corrections to this simple exponential decay 
become important, when F has the same order of magnitude as the mean level spacing 
of the background Hamiltonian. 

In Ref. [6] corrections to the FGR-law, were calculated in a mesoscopic system, 
where the special state is sitting on a quantum dot which is weakly coupled to a reservoir. 
Due to these corrections the decay of the electronic state in the dot is never complete. 
Instead the system preserves a memory of the original state. The corrections to FGR, 
which are in nature very similar to the weak localisation correction to classical transport, 
give rise to non-ergodic behavior. In the same spirit, in Ref. [7] weak localisation 
corrections to the FGR behavior were calculated using new semiclassical techniques, 
which became available recently 0,19]. 

Corrections to the FGR-law can be addressed in a more generic setting within a 
random matrix model. In this model a special state = \s) couples to a large 

reservoir of states via randomly chosen coupling parameters V^, fi = 1, . . . N. The 
coupling parameters are chosen either real or complex. The dynamics of the reservoir 
is either chaotic, modelled by a random matrix chosen from a Gaussian RMT ensemble 
or regular, with Poissonian eigenvalue statistics. This RMT model has been addressed 
in plU El] and corrections to FGR regime were found. 

Expanding the special state in eigenstates of the Hamiltonian survival probability 
can be written as 



P(t) 



w)\m)\ 2 - 



(i) 




(2) 
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For uncorrelated eigenvalues and expansion coefficients (s\m) the ensemble average 
(denoted by a bar) can be performed for both sums separately. As a result 

W) = \W)\\ (3) 
where p(t) is the Fourier transform of the local density of states LDOS 

p(E) = J2\(s\n)\ 2 5(E-E n + E s ) (4) 

n 

around the energy E s of the special state. As was pointed out already by Weisskopf and 
Wigner [12J the smooth part of p(E) is under very general assumptions of Lorentzian 
shape 

i r/2 

n(E-E s y + (T/2f 

from which FGR is recovered. We call the approximation implied in Eq. D(rude) 
B(oltzmann)- approximation, due to its formal similarity to the approximations made in 
the derivation of the Drude Boltzmann law of conductivity [5]. 

Assuming a constant mean level spacing D in an energy region around the special 
state, from Eq. (j2J) one obtains P(oo) = IPR, where 

IPR = ^|(s|m)| 4 = D I dEp 2 (E) (6) 

m 

is the inverse participation ratio of the special state in the basis of the eigenstates of 
the full Hamiltonian. In the DB-approximation p 2 ~ p 2 and the mean IPR can be 
estimated by IPR ~ D/(tvT), which shows that survival probability will not decay to 
zero, if T and D are of the same order of magnitude. 

Following Ref. [10] one can obtain corrections to the DB approximation due to 
energy correlations as follows: Writing the expansion coefficients |(s|n)| 2 = p(E n ) + 
5p(E n ) as the sum of a smooth function of E n and of a fluctuating part, the averaged 
survival probability can be written as the sum of two contributions as well 

Pjt) = J dEdE'p(E) p(E')e l{E - E ' )t R 2 {{E - E')/D) + ~5P(f) , (7) 

where R 2 = m S(E — E n )6(E' — E m ) is the averaged energy-energy correlator of the 
background Hamiltonian. The approximation made in Ref. [10] consists in neglecting 



the fluctuating part SP(t). Using standard results of random matrix theory [13] R 2 is 
the sum of a 5 like , a connected and of an unconnected contribution. Likewise, survival 
probability is written as a sum of three contributions 

oo 

Tit) ~ e~ Tt + ^-- [ dt'e-^-^b^t') (8) 
7rr J 

— oo 

where b 2 is the two-level form factor. The last term in Eq. (jSJ) accounts for the energy 
correlations of the bath. The result (jSJ) gives an intuitive insight how energy correlations 
of the background Hamiltonian give rise to corrections of the FGR-law and as we will 
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see for strong coupling T ^> D and for a chaotic background it predicts qualitatively 
the correct behavior. 

However, it is easy to see that for small T or in the case where correlations are 
absent Eq. ([8]) is not correct even qualitatively (for instance the saturation value P(oo) 
= D/iirY) exceeds one for Y < D/ti). 

Thus the question, whether Eq. ([8]) describes sufficiently the weak localisation 
corrections to the DB-approximation has to be answered negatively. In the present 
work we therefore calculate P(t) for the random matrix model mentioned above exactly 
for a chaotic as well as for a regular background. We will see that the exact result 
differs qualitatively from the predictions of Eq. ([8]). For instance, we find that a revival 
of survival probability, predicted by Eq. ([8]) only for fairly strong couplings and for a 
chaotic background, occurs for weak coupling and for a regular background as well. 
More general, the energy statistics of the background turn out to have little influence 
on the survival probability. Instead, we find that the nature of the coupling coefficient 
is crucial. As a rule of thumb, for constant mean coupling strength survival probability 
is always lower for complex coupling coefficients than for real ones. 

We provide exact analytic expressions for the average IPR, which interpolate the 
power law decay for strong couplings to the small coupling regime. 

On the technical level, we use powerful results for averages over characteristic 
polynomials, which have become available recently [14] . This allows us to circumvent 
a long and complicated supersymmetric calculation. This elegant shortcut is possible 
for complex couplings of the doorway state to the background, where we derive exact 
analytic results. For real coupling coefficients we resort to numerics. 

The paper is organized as follows: In Sec. [2] and in Sec. [3] we define the random 
matrix model and fix the notation. In Sec. H] we outline how the Lorentz shape of the 
LDOS comes about in the present random matrix model. The calculation of survival 
probability for regular background as well as for a GUE and for a GOE background is 
presented in Sec. 15.21 Finally the results are discussed and summarized in Sec. El 

2. Definition of the Doorway model 

The model to be discussed here stems from nuclear physics [15J and is also often used 
in other fields jT6| . For the convenience of the reader and to define our notation, we 
compile its salient features. 

The total Hamiltonian H consists of three parts, the Hamiltonian H s for the 
doorway states, the Hamiltonian describing the N background states, where N will 
eventually be taken to infinity, and the interaction V coupling the two classes of states. 
Often, there is only one relevant doorway state or the spacing between the doorway 
states is much larger than their spreading widths. In the present work we focus on this 
situation, leaving the interesting case of many doorway states to future work. Hence, 
we have 



H = H S + H b + V 
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H = H + V 

N N 

= E.\s)(s\ + ^2E v \b v )(b v \ +^(v v \s)(bv\+h.c.} . (9) 

For the matrix elements of the interaction, we make the assumptions (b v \V\b^) = and 
(6 v |V|s) = V v for any fi, v. 

Resembling the situation in most systems, we put the doorway state \s) in the center 
of the background spectrum. It interacts with the N surrounding states. Without loss 
of generality, we may set E s = 0. The eigenequations for the uncoupled Hamiltonian 
Hq are then 

H s \s) = and H b \b v ) = E u \b u ) . (10) 

We assume that the interaction matrix elements are Gaussian distributed random 
variables. We distinguish the two cases of complex (/3 = 2) or real (j3 = 1) matrix 
elements V v . Introducing the iV-component vector V, the corresponding distribution is 

As discussed in the introduction, we are interested in the situation where the mean 
coupling strength is of the same order of magnitude as the mean level spacing. We 
define the dimensionless parameter 



_ yW) _ v 

x ~ ~7Wb~ - D ' (12) 

where D is the mean level spacing of the background states in the center of the 
band [151 [16]. The distribution Pi(V) is chosen such that A is independent of j3. 

We distinguish two cases for the background dynamics. Regular dynamics of the 
background Hamiltonian Hf, is modelled by uncorrelated eigenvalues 

N 

P b (H b ) = Y[p b {E v ) . (13) 

u=l 

Chaotic dynamics is modelled by Gaussian random matrix ensembles given by the 
distribution function 

JV ,3 b N(N-l) 
A> V fPb\ 4 ( A>._ TT 2 



where H b is either real symmetric {(3 b = 1) modelling time reversal invariant background 
dynamics, or Hermitean (f3 b = 2), modelling background dynamics with broken time 
reversal invariance. The probability distribution ( fl"4"l) yields a mean level spacing 
D = a/ tt 2 /2N in the band center, which is independent of f3 b . We denote the average 
over both the interaction matrix elements and the background Hamiltonian by a bar 

(TO = j d[H b }P b (H b ) J d\y\Pi(V)(...) . (15) 
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3. Fidelity and Survival probability of the Doorway state 

We define the echo operator 

M\ = e iHt e- lHot . (16) 

Fidelity amplitude f\(t) is defined as the expectation value of the echo operator with 
respect to a given initial state. Here we are interested in the Doorway state \s) as initial 
state, i.e. an eigenstate of the unperturbed system. Since E s \s) =0 the average fidelity 
amplitude can be written as 

Mt) = (s\M x (t)\s) = (s\e^ m \s). (17) 

As mentioned in the introduction, fidelity amplitude is then the Fourier transform of 
the local density of states p(E). Likewise, fidelity (often called Loschmidt echo) F\(t), 
defined as the modulus square of the fidelity amplitude, becomes identical with the 
survival probability P(t) 

F x (t) = (s\e- iHt \s) (s\e iHt \s} = P(t) . (18) 

In the following we mainly stick with the notion of survival probability, keeping in mind 
that in the present situation fidelity and survival probability are synonyms. 

When we expand fidelity in eigenstates of the full Hamiltonian according to Eq. (j2J), 
we see that in the limit of infinite large times the fidelity approaches the inverse 
participation ratio IPR of the special state in the basis of the eigenstates of the system. 
We therefore also define the mean inverse participation ratio 

IPRI = ^|(s|m>| 4 = F\^o), (19) 

rn 

where the sum goes over exact eigenstates of the full Hamiltonian H. The task is 
to calculate f\(t) and F\(t) exactly in the large N limit for various choices of the 
background Hamiltonian and in particular the corrections to the DB-approximation. 



4. Calculation of the mean local density of states 



We first prove that the average local density of states p(E) in the present case indeed 
has the Lorentz shape. We write 

p(E) = ^|( S |m)| 2 5(E - E m ) 



-Im 



7T 



7T det(# 



H-E-ie 
det(H - E) 



E-ie)' (20) 
where we used Kramer's rule in the step from the second to the third equation. For a 
chaotic background the average can be taken most conveniently by a mapping onto a 
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super symmetric matrix model. Using standard steps [T7] one arrives at 



p(E) = hm J d[a) exp (-^ Str(a + Ef^j SdeU 



2 (*) 



det"^ 2 



+ E 



0"BB 



(21) 



where the bar now denotes an average over the coupling coefficients only. In Eq. (l2Tj) o 
is a 2 x 2 (GUE, (3b = 2) respectively 4x4 (GOE, (3 b = 1) supermatrix of the form 

, GUE 




GOE . 



(22) 



The matrix entries in Latin letters denote real commuting integration variables. The 
matrix entries in Greek letters denote complex anticommuting integration variables. The 
infinitesimal volume elements d[a] are products of the differentials of all independent 
integration variables. The integration domain of the real commuting variables is the 
real axis. The so-called Boson-Boson block cxbb is the upper left block of commuting 
variables. The matrix integral Eq. (|2T|) can be solved in one step with a saddlepoint 
approximation. For energies close to the center of the band the saddle points are 
a ~ 



±i^N/2 and thus 



p{E) 



—lm.\E - 

71 V 



(23) 



This result shows that LDOS has the form of a 5-spike unless the perturbation is 
classically small, i. e. of the order of the mean level spacing. The Gaussian average over 
the coupling coefficients finally yields the Lorentz distribution 

i r/2 



P(E) 



7iE 2 + (T/2) 2 



with the crucial relation 
T = 2rr\ 2 D 



(24) 



(25) 



between spreading width, mean perturbation strength and mean level spacing. For a 
regular environment the LDOS was calculated for instance in [18] yielding the same 
result. 



5. Calculation of mean Fidelity /Survival Probability 

We now turn to the main task: the calculation of survival probability of the doorway 
state. In order to calculate the mean fidelity /survival probability, we write F\(t) as 

= ^2 / dEl j dE * ^ - P( E * + ie )P( E * ~ lt ) 
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dEi I dE 2 exp(i(E 1 - E 2 )t) 

det(H -E 1 ) det(H -E 2 ) 
det(H - E 1 -ie)det(H - E 2 + ie) ' 1 j 

where we used again Kramer's rule. Evaluation of the determinant in the denominator 
yields 



F x (t) = J dE x J dE 2 exp (i(E 1 - E 2 )t) 8 ^E ± + 



\V U 



Ea — Ei 



S[E 2 + J2 



E,, — Ei 



(27) 



We observe that F is normalized by F\(0) = 1 (see [TQl [19] ) . This allows us to extract 
at this point of the calculation the constant term F\(0) from the integral and to average 
over F\(t) — 1 instead of F(t) directly. After a Fourier transformation of the delta 
distributions we find for the average 

Wt) = l + JdE 1 JdE 2 J^-J ^[exp^-i^))-!] 



exp 



-i k -^{E l -E 2 )-r k -±±^{E l + E 2 ] 



exp —i 



^ TP - 



V n 2 ko 



(28) 



E u — E-\ — 4 E,. — En 

Since on the unfolded scale the mean level spacing of the background is constant, we 
can assume that the unfolded average does not depend on E\ + E 2 . This allows us to 
perform the integral over the mean energy (E\ + E 2 )/2 and, trivially, over {k\ + k 2 ). We 
find 



/f dk 
ds I — [ ex P ( iTS )) ~ !] ex P (-iks) R(k, s) 



where 



R(k, s) = exp j — i 



Kl'ks 



El - (Ds/2y 



(29) 



(30) 



We introduced the dimensionless time r = Dt, measured in units of Heisenberg time 
7h = D~ x . It is useful to take the average over the Gaussian distributed coupling 
coefficients at this stage of the calculation 



R(k,s) 



x ' d[V,f 



n 



\ \ (2nv 2 / p)?/ 2 



exp 



ft 



2v 2 



1 + 



2iksD 2 \ 2 



(3{E 2 - {Ds/2f 



det(H 2 - (Ds/2)< 



det(H 2 - (Ds/2) 2 + iksD 2 X 2 



p/2 



(31) 
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where the bar denotes the average over the background Hamiltonian Hi, only, which 
has still to be performed. In the expression ([3T1) it becomes evident why the case of 
real coupling coefficients G R, (3 = 1 is analytically more difficult than the case of 
complex coupling (3 = 2. For (3 = 2 the expression (I3T1) is an average over a rational ratio 
of products of characteristic polynomials. Much information has been gathered about 
these averages in the last decades [201 1211 (221 [231 1211 [251 [26] , whereas little is known 
about averages over irrational functions of characteristic polynomials as encountered in 
the case (3 = 1. 

In the subsequent analysis we therefore restrict ourselves to complex coupling and 
set j3 = 2 from now on. In the case of real coupling we recur to numerics. We distinguish 
the cases of a regular background -Rp isson(^, s), and GOE or GUE distributed chaotic 
background RGOE{k,s), -Rgue(&, s). We are able to calculate all three averages exactly 
in the large N limit. 

5.1. Survival Probability for a regular Background 

We first consider a regular, Poisson distributed, background. For Poisson distributed 
eigenvalues the average in Eq. ( 13T1) becomes a product 

#Poisson(M =r(k,s) N (32) 

r / x 2 -(s/2) 2 \ f3/2 

r(k,s) =Dj dxPb (D X ) ( l2 _ (s/2) ^ faAV3 ) • (33) 

The universal final result should be independent of the distribution of the eigenvalues 
of the background Hamiltonian. The simplest choice for this distribution is 





\x\ 


{o, 


\x\ 



pb( E ) = 77^' Llr.V/o (34) 



where D = 1/ y/N and \/N = ND is the length of the background spectrum. 

For (3 = 2 the integral (BSD can be evaluated in the large N limit by the residue 
theorem 

r(k,s) =1 . 1 " 1 +0\m) , 35 

/ 7 \ ( 2ttA 2 |A;||s| \ . . 

^Poisson (k,s) = exp . (36) 

V y s z — Aiks\ z J 

Using this result together with Eq. ( 1321 and Eq. ( 1291) we find 

oo 

^) - 1 = Re y — [cos(rs) - 1] 



o 



which is almost our final result. However the remaining double integral is numerically 
difficult due to the oscillatory terms. We proceed by rotating the contour of k integration 
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on the negative imaginary axis. Introducing the new integration variable x = ik/s one 
arrives at 

oo 1/4A 2 

— 7 — r 1 f if, I 2\ ( 27rA 2 xs \ 

F\(t) = 1 H — / ass cos(rs) — 1 / ax exp — xs ) sin = 

7T J 1 ^ J V ^ W1-AX\ 2 J 



(38) 



The s integration can now be performed without difficulties. The final result is 



i 

A f dx _*ze£a!_ f 7T / t 2 x 2 , ( 7lX 2 T 



F x (t) = 1 H = / —=e «u-«o ^ - e s~ cosh - 1 



o 



2r t 2 a 2 , / 7rA 2 r . . 
e — smh )■ . (39) 



x \1 — x ^ 

The remaining integral does not seem to have a simple analytic solution. As expected 



IPRa is not zero but saturates at a finite value, given by 



7r 3 A /7rA\ 2 „ //tA 



IPR A = l-^-exp^ T j erfc K-J • (4°) 
This function behaves for small/large values of A as follows 
f 1 - tt 3 / 2 A for A«l 

In Figure [1] survival probability is plotted on a logarithmic and on a linear scale as 
a function of time in units of Heisenberg time for three different values of the mean 
coupling strength A = 0.1, 0.5 and 1 corresponding to a spreading width Y/D ~ 0.06, 1.5 
and 6.3. It is seen that the survival probability reaches a minimum and increases 
afterwards to its saturation value given in Eq. (HU|) . The time evolution of survival 
probability splits into three regimes: for t tr/Y fidelity follows the FGR law, for 
t 3> Tn/r survival probability has approached its saturation value and is approximately 
constant, in the region t ~ th/T the time evolution is a complicated smooth function, 
which interpolates between the two limiting regimes. 

In Figure [1] also the curves obtained from Eq. (jSJ) are plotted. It is seen that for a 
Poisson distributed spectrum of the background, Eq. (jSJ) is a rather poor approximation 
of the exact curve. In particular it predicts no revival of survival probability and the 
saturation value is underestimated by a factor four for large coupling strength A. 

5.2. Survival probability for a chaotic Background: GUE 

Using the formulae of Theorem 1.3.2 of Ref. [U] we find for R(k, s) for a chaotic 
background Hamiltonian with broken time reversal invariance 

-Rgue(&, s) = exp (— iirsgn (ks)V s 2 — AiksX 2 ^ 

\ For the convenience of the reader we provide Theorem 1.3.2 of Ref. [14] in the form needed here in 
|Appendix A| 
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Figure 1. Plot of Eq. (39]) for the values A = 0.1 (red thick line), A = 0.5 (blue thick 
line) and A = 1 (green thick line) on a logarithmic scale (right) and on a linear scale 
(left). The curves obtained by Fermi's golden rule are depicted by thinner dotted lines 
in all three cases. For A = 1 and for A = 0.5 also the curves obtained from Eq. ([8]) are 
plotted (thin dashed lines). 



cos(7rs) + isgn (ks) sin(7rs) 



s - 2ik\ 2 
y/s 2 - AiksX 2 



(42) 



We use that R(—k,s) = R*(k,s) and rotate the contour of the fc-integral as in the 
Poisson case to the negative (s > 0) or positive (s < 0) imaginary axis. We find for the 
averaged survival probability 

oo 1/4A 2 

F\(t) = H — J dss [cos(ts) — 1] J cfeexp(— xs 2 ) 



sin ( ns^/l — AxX 2 ) cos(7rs 





COS I 71 



syl — AxX 2 ) sin(7rs 



2X 2 x 



VI - 4xA 2 



(43) 

The s integration can be performed in a tedious but straightforward way The final 
result is again an integral expression 



F x (r) = 1 + 



A 



dx 



2v^ Jo s/x\/l^x 

\ 2 T 2 

II - ' ' 

X 



exp 



7T 2 X 2 Wl 



X 



X 



+t exp 
A 



1 — exp 
X 2 r 2 



x 



sinh 



cosh 
2ttX 2 tW. 



2ttX 2 tW+ 



x 



X 



dx 



X 2 r 2 



exp 



7lW- 

+t exp 



1 — exp 

x 



sinh 



cosh 



2ttX 2 tW_ 



x 



7T 2 X 2 W 2 

X 

2ttX 2 tW. 



X 



(44) 
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2 4 6 8 10 5 10 15 20 

Figure 2. Plot of Eq. ([44]) for the values A = 0.1 (red thick line), A = 0.5 (blue thick 
line) and A = 1 (green thick line). The curves obtained by Fermi's golden rule are 
depicted by thinner dashed lines in all three cases. For A = 1 and for A = 0.5 also the 
curves obtained from Eq. |(5J) are plotted (thin dashed lines). 



where 

W± = l±VT^x. (45) 

In Figure [2] survival probability is plotted on a logarithmic and on a linear scale as 
a function of time in units of Heisenberg time for three different values of the mean 
coupling strength A = 0.1, 0.5 and 1. These values of A correspond to a spreading width 
T/D « 0.06, 1.5 and 6.3. It is seen that qualitatively the curves are quite similar to the 
ones obtained for a regular environment. For a GUE background energy correlations are 
present and a revival of survival probability is predicted by the estimation (jSJ). However 
we notice that a revival occurs also for small values of coupling strength like A = 0.1. 

The averaged inverse participation ratio is easily obtained from Eq. (jUj) by taking 
the limit r — > oo 

Av^F f dx ( tt 2 A 2 (2-x) 



IPR A = 1 — / — . exp 



2 J y/xy/l - X V x 



cosh I I + v 1 — x sinh I I . (46) 

The remaining integral can be simplified and expressed in terms of complementary error 
functions akin to Eq. (j401 



IPRa = 1 - 



, 00 9 

2vr 4 A 4 f e~ u du 



^ J (M 2 + 7T 2 A 2 ) 2 







7r 2 A 2 exp ^2 A 2J erfc(7rA ) . ( 47 ) 



In the limits of small and large A we obtain 

1 - tt 3 / 2 A for A«l 



1 (48) 
-r-— for A > 1 . y 1 

77 A 
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The asymptotic value of IPRa is half the value obtained for a regular environment but 
twice the value predicted by the Boltzmann-Drude approximation. 

5.3. Survival Probability for a chaotic Background: GOE 



Using Theorem 1.3.1 of Ref . . for details see | Appendix A[ an expression for -Rgoe(^, s) 
can be derived. We find that 

-RgoeO, s) = -RgueO, s) + Rzdd(k, s) 

„ . — iAsgn (ks)k 2 sX 4: ( d sin(7rs)\ 

^(fc, s) = — f=r- i x 

Vs 2 - AtksX 2 \ds s J 



oo 



exp ^— z'7rsgn (ks)Vs 2 — AiksX 2 tj — . (49) 



t 



i 



In the same fashion survival probability splits into two parts 



7 A,GOE — -F\,GUE + -^A.add 



A, add 



/r (ik 
ds I — [exp {its)) - 1] exp (-iks) R ad d(k, s) . (50) 

The additional contribution to the survival probability F\ A dd plays the role of a 
Cooperon contribution. In a tedious but straightforward calculation we find for -FA >a dd 
an expression as a double integral 

1 (X) 



A, add 



jdxj dt g ^— (H(t) + H(-r) - 217(0)) 



o 1 
H(r) = e^l^M 2 ) 



4A 2 tt 1\ . [2tx X 2 W{t) 



+ - W(t) sinh 



X IT J \ X 

(vr 2 + W(t) 2 ) cosh 

x \ x 

W(t) =r + ntVT^ , (51) 

which can be evaluated numerically without problems. Likewise the average of the IPR 
obtains an additional contribution 



IPRa, GOE — IPRa,GUE + IPRA,add 
1 oo 

IPR^ = - f dx [ dt ^^B . (52) 
J J At a/1 - x 

o 1 

On the left hand side of Fig. [3] survival probability as given by Eq. (150]) is plotted for 
different coupling strength A = 0.1, 0.2 and 0.5 (full lines). A comparison with the 
corresponding curves for a GUE background (dotted lines) shows that the difference is 
minimal. Whether or not the background dynamics is time reversal invariant or not 
has no influence of the decay of the special state. Nevertheless, it is interesting to look 
at the Cooperon contribution F\ ia dd separately. It is plotted on the right hand side of 
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Figure 3. Left: Plot of Eq. ([50]) and Eq. JSTJ for the values A = 0.1 (full black line), 
A = 0.2 ( full red line) and A = 0.5 ( full green line). For comparison the curves for 
the GUE are plotted for the same values of A with dotted lines. Right: Plot of the 
additional "Cooperon" contribution to survival probability according to Eq. (|5Tj) for 
the three values A = 0.1 (black line), A = 0.2 (red line) and A = 0.5 (green line). 

Fig. [3] for the same values as before. We see that the contribution is small compared to 
F\,gve- Surprisingly, we see that it is oscillating a few times with a frequency oc 1/A 
before reaching its saturation value IPR^add which is a non-monotonous function of A. 
This contribution to the total mean IPR is depicted in the inset of Fig. H] as a function 
of the coupling strength A. 

5.4- Comparison 

We were able to calculate time evolution of survival probability for a complex coupling 
of the prepared state to a Poisson, GUE or GOE environment. As explained before a 
similar calculation is by now not possible for real coupling coefficents. In the latter case 
we resort to numerics. 

In Fig. H] time evolution of survival probability is plotted for regular and GUE 
background dynamics for real and for complex coupling coefficients (the difference 
between a GOE background and a GUE background is almost invisible on the scale 
used for the plots). We see that increased background complexity reduces overall 
survival probability. This is in agreement with standard perturbative arguments [21 E]. 
The difference between real and complex coupling coefficients is of the same order of 
magnitude as the difference between regular and chaotic background dynamics. This is 
surprising inasmuch time reversal symmetry breaking in the background has practically 
no influence on survival probability. 

In Figj5]the average inverse participation ratio is plotted for a complex coupling to 
a Poissonian background (as given by Eq. (I40p ) and to a GUE background (as given by 
Eq. (H7J) as a function of coupling strength A. We see that for A > 1 the mean IPR is well 
approximated by its asymptotic form, Eq. (j4"T|) respectively Eq. (|4"8l) . In the inset the 
additional contribution for a GOE background IPRA, a dd is plotted. Although negligible 
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Figure 4. Left: Comparison of survival probability for real coupling to a Poissonian 
background (full red line), to a GUE background (full black line), for complex coupling 
to a Poissonian background (dashed red line) and to a GUE background (dashed black 
line). The coupling strength is A — 0.5. The inset shows the same quantities for 
coupling strength A = 0.1. 



IPR x 




Figure 5. Left: Plot of the average inverse participation ratio IPRa as a function of 
the dimensionless mean coupling strength A for complex interaction with a regular 
background (green curve), Eq. (|40j) . and for complex interaction with a GUE 
background (red curve), Eq. (j47|) . The blue dashed lines show the asymptotic behavior 
for large A. The inset shows the additional contribution for a GOE background 
IPR^add as a function of A. 



for practical purposes, it is interesting to see that this contribution is a non-monotonous 
function of A. It vanishes for A = and for A = oo. It obtains its maximum for A ~ 0.5. 
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6. Discussion and Summary 

We calculated exactly survival probability and fidelity amplitude for a special state, 
which is weakly coupled to a random matrix environment. Whereas fidelity amplitude 
decays exponentially according to Fermi's golden rule, survival probability shows a rich 
behavior. We found a revival of survival probability after a characteristic time which 
increases with decreasing coupling strength and a saturation of survival probability at a 
value given by the mean IPR of the special state in the basis of the interacting system. 
Our exact results largely improve existing estimates [10] of these quantities. 

We were able to derive analytically the full A-dependence of the IPR in the small 
coupling regime, where the approximation IPR oc A -2 becomes bad. It turned out 
that even in the strong coupling limit IPR is largely underestimated by the Drude- 
Boltzmann approximation (by a factor four for a regular environment and by a factor 
two for a chaotic environment). 

Revival of survival probability was found for all types of background complexity. 
The fact that it occurs also for a regular background encumbers an explanation of the 
revival by spectral correlations of the background energy-levels as put forward in [10J. 
The revival is quite different in nature to the fidelity revivals reported in [271 1281 29J . 
There a global perturbation and fidelity of a random state was considered and a revival 
of fidelity at Heisenberg time was found. An explanation of this phenomenon was 
given by the rigidity of the spectrum of a background with chaotic dynamics. Such 
an explanation obviously fails in the present case, since revival occurs even in the 
absence of energy correlations of the background. The behavior rather resembles the 
overdamped oscillations in a two-level system, which is coupled to a non-Markovian 
heat-bath (for instance one or more spin baths [301 [31]) • A possible explanation is that 
for small couplings the system effectively reduces to a two-level system involving only 
the Doorway state and its nearest neighbor (in energy). This two-level system itself is 
then strongly coupled to the remaining background states. 

Survival probability is susceptible against changes in the background dynamics from 
regularity to chaoticity, this is in agreement with the original arguments of Peres [32J. 
It is not sensitive against time reversal symmetry breaking in a chaotic background 
dynamics. 

Whereas fidelity amplitude has been calculated exactly in various random matrix 
models, this is the first exact calculation of Fidelity in a random matrix setting. This 
was possible, due to advances in the calculation of ensemble averages of characteristic 
polynomials in the last years [2QIEI1E21E31E21EE1III1EE]. The exact results are limited 
to the case of a complex (time reversal invariance breaking) coupling to the background. 
A similar calculation for real (time reversal invariant) coupling would require knowledge 
of the averages of non-rational functions of characteristic polynomials. 
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Appendix A. Theorems 1.3.1 and 1.3.2 of Ref. [14] 

For a chaotic background the ensemble average can be done using a result by Borodin 
and Strahov [TJ]. We restate Theorem 1.3.2/1.3.1 of Ref. Pj concerning the GUE/GOE 
ensemble average of a ratio of an arbitrary number of characteristic polynomials in a 
form adapted to our purposes. 

We first state Theorem 1.3.2 concerning the GUE: Define the multivariate function 
C(a,/3) of 2n + 2m variables a^, a^, 8f, 1 < k < n, 1 < I < m with a^,j3[ G C 
and <y^,/3^ G C\M as the ratio of an arbitrary number of characteristic polynomials 

and the average over N x N random matrices chosen from GUE 

d[H]P b {...) , (A.2) 



where the distribution P b is given by Eq. (jT4"|) . with (3 = 2. Moreover, define 
7 = (n + m) 2 + (n — m) and the Vandermonde Determinant A n (a) = Hi — a j)- 
Then the following identity holds 



lim C(a/v / 2iV,/?/v / 2iV) 



v ' ,r/ ' v ' A n (a-)A n (a+)A ro (/3-)A m (/3+) 

det[S (2 \a~,fJ + \fJ-,a + )} , (A.3) 

where S^(a~ , , a + ) is a n + m matrix with rows parametrized by elements a^, 

and columns parametrized by elements (3^, a\ and with matrix elements [§] 

m\ 1 sin(a~ — 8~) 

7T a p - 8 q 

( _ exp i ( Q j--a p -) j + q 

S( 2 )( a ;,«+)= ^ 9 (A.5) 



expi(/3+-/3~) 



Ima+ < 0, 



Im/3+ > 0, 



§ In the second line of Eq (|A.7j) there is a minus sign changed as compared to the original theorem of 
Ref. [14] . which is apparently wrong. 
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Pp -a q 



. Im/3+ > 0, Ima+ < 0, 



fi^0?,a+)=27ri^ e ^f_~f\ Im/3+ < 0, Ima+ > 0, (A.7) 

0, in all other cases. 

Since the mean level spacing is given by D = n/\ / 2N, we find that R{k, s) is exactly of 
the form (1A.3j) with n = m = 1 and with 

TVS ITS 



^Vs 2 - 4zA;sA 2 , /?+ = ~Vs 2 - 4iks\ 2 . (A.8) 
This yields Eq. (p|). 

We now turn to Theorem 1.3.1 concerning the GOE: Define the multivariate 
function C(a,(3) of n + m variables a^, (3i, 1 < k < n, 1 < I < m with a& G C 
and $ G C\R as the ratio of an arbitrary number of characteristic polynomials 

and the average over 2iV x 2N random matrices chosen from GOE 



a 



(...)=/ 4tf]A(-..) , (A.io) 



where the distribution P^ is given by Eq. ffl4|) . with (3 = 1. Then the following identity 
holdil 

lim C(a/v / 4iV,/?/v / 4iV) = E% ~ 

Pl[S {1 \a,(3\a,(3)] , (A.ll) 

where S^(a~, (3 + \(3~ , a + ) is a skew-symmetric n + m matrix with rows and columns 
parametrized by elements a, /3 and with matrix elements 

Q(i)t \ 1 9 sin(a p - q g ) 

S i; (ap,a 9 )= , (A.12) 

7r octi a p — a q 



expi(f3 q -a p ) 



lm(3 q > 0, 



Setting 



s" M )={ ^s^( Im ,; <0 ; (A.i3) 

jH-oo exp^-ft)*) ^ Im/3p > Q) < 0> 

SV((3 p ,P q ) =2m{ _j+^ ESm=m dt, TmP p < 0, Im& > 0, (A.14) 

0, in all other cases. 



ITS ITS 

1 2 2 



ft = ^Vs 2 - AiksX 2 , ft = -^Vs 2 -AiksX 2 . (A.15) 
in Eq. flATTTj) yields Eq. (|49|) . 

|| a scaling factor y2 seems to be wrong in Ref. [H] 
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